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Reissner-Nordstrom black holes have two static regions: r > ro and < r < r^, where 
and To are the inner and outer horizon radii. The stability of the exterior static region has been 
established long time ago. In this work we prove that the interior static region is unstable under 
linear gravitational perturbations, by showing that field perturbations compactly supported within 
this region will generically excite a mode that grows exponentially in time. This result gives an 
alternative reason to mass inflation to consider the space time extension beyond the Cauchy horizon 
. as physically irrelevant, and thus provides support to the strong cosmic censorship conjecture, which 

j—i ' is also backed by recent evidence of a linear gravitational instability in the interior region of Kerr 

I black holes found by the authors. The use of intertwiners to solve for the evolution of initial data 

■ plays a key role, and adapts without change to the case of super-extremal Reissner-Nordstrom black 
I ' holes, allowing to complete the proof of the linear instability of this naked singularity. A particular 

^ ' intertwiner is found such that the intertwined Zerilli field has a geometrical meaning -it is the first 

, order variation of a particular Riemann tensor invariant-. Using this, calculations can be carried 

' out explicitely for every harmonic number. 

I— l' PACS numbers: 04.50.-|-h,04.20.-q,04.70.-s, 04.30.-w 

o ; 
cr. 

\-i ! I. INTRODUCTION 

In the course of a program [1-4] to study the stability under linear gravitational perturbations of the most notable 
nakedly singular solutions of Einstein's field equations, namely, negative mass Schwarzschild's solution [1, 4, 6], 
jft , IQI > M > Reissner-Nordstrom space-time [2] and \J\ > Kerr space-time [2, 3] (see also [7]), we noticed that 
' the stationary interior region beyond the inner horizon of a Kerr black hole is unstable [3, 5]. The existence of an 
I , initially bounded and exponentially growing solution of Teukolsky equations in the "super-extreme" case \J\ > M^, of 

■ which some numerical evidence had been given earlier in [2] , was established in [3] , where it was shown that there are 
' I actually infinitely many unstable modes. It was also shown in [3], that the stationary region beyond the inner horizon 

• of a Kerr black hole (i.e. \J\ < M"^) is linearly unstable under gravitational perturbations. These results show that 
linear perturbation theory is a valuable tool to study not only weak cosmic censorship (impossibility of formation of 
T-H ■ naked singularities), but also strong cosmic censorship (impossibility of formation of Cauchy horizons). 
ILjl ' For the Kerr spacetime, an explicit expression for the unstable modes is not given in [3] , since they involve solutions 
of complicated second order ordinary differential equations, which can at best be written in terms of Heun functions, 
providing little extra information. This, added to the complexity of the reconstruction of the perturbed metric from 
5^ , a solution of Teukolski's equations, makes it extremely difficult to evaluate the physical meaningfulness of these 
unstable modes. 
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The situation is different for the negative mass Schwarzschild and the super-extremal Reissner-Nordstrom space- 
times, where explicit expressions for some unstable modes, which involve only elementary functions, are given in [1, 4] 
and [2] respectively. Since the metric reconstruction process in the spherically symmetric case is much simpler, it is 
possible to study the effect of perturbations on the singularity, by calculating the perturbed Riemann tensor invariants. 
We use these to select appropriate boundary conditions at the singularity that ensure the self-consistency of the linear 
perturbation scheme, by requiring that curvature scalars do not get corrections that diverge faster at the singularity 
than the zeroth order term. As an example, in the case of the Schwarzschild negative mass naked singularity there are 
infinitely many possible boundary conditions at the singularity, parameterized by [1, 6], only one of which satisfies 
the above requirement. Thus, besides assuring the self-consistency of the perturbative treatment, the above procedure 
solves the problem of having a unique, well defined evolution of perturbations in a non globally hyperbolic background. 

The unstable modes in [1] were recognized by Cardoso and Cavaglia [7] to correspond to Chandrasekar's "algebraic 
special" (AS) solutions of the linearized Einstein's equations [9, 10]. This observation hinted in the right direction 
where to look for unstable modes of super-extremal Reissner-Nordstrom black holes [2, 7]. As shown in [2], some of 
the Reissner-Nordstrom AS modes grow exponentially in time while keeping appropriate spatial boundary conditions 
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in the super-extremal case, as happens in the negative mass Schwarzschild case. With no exception, the AS sohitions 
are irrelevant to the stability problem of the exterior region of black holes, since they do not satisfy suitable boundary 
conditions, a probable reason why they remained unnoticed for such a long time. For the Kerr solution, the AS 
modes do not satisfy appropriate boundary conditions, neither for the black hole stationary regions, nor for the 
naked singularity. However, it was proved in [3] that unstable modes exist for every harmonic (i.e., spin weighted 
spheroidal harmonic) of the Teukolski equations in the supc!r-c!xtrcmc case. Moreover, in [3] a connection was 
spot between the unstable modes of Kerr naked singularities, and unstable modes for the interior stationary region 
(r < rj) of a Kerr black hole. Given the similarities in the structures of the maximal analytic extensions of Kerr and 
Reissner-Nordstrom black holes, and the fact that both are affected by Cauchy horizon issues, one is naturally led to 
ask whether the interior region of a Reissner-Nordstrom black hole is also unstable. In this paper we show that this 
is the case. We give a detailed proof of the instability of the inner region of a Reissner-Nordstrom black hole under 
linear gravitational perturbations initially restricted to a compact subset of the inner region. 

This provides an alternative reason to the mass inflation mechanism to disregard the extension of the space time 
manifold beyond the Cauchy horizon: since the inner region is static but unstable, it cannot be the endpoint of an 
evolving space time. A similar result for the Kerr black hole would imply cutting off the inner region of this space 
time, which has closed time like curves and other pathologies. 

We will concentrate on type one (also called "gravitational" , as opposed to typo two or "electromagnetic" ) scalar 
(also called "polar") linear perturbations of the metric and electromagnetic fields around the Reissner-Nordstrom 
solution, since it is in this sector that we have found explicit unstable modes. The linearized Einstein's equations for 
these type of perturbations can be reduced to a 1 + 1 wave equation on a field in a semi-inflnite domain bounded 
by the singularity world- line, with a time independent potential (Zerilli's equation [11-13]). This formalism was used 
to prove the stability of the exterior static region [11] of the Reissner-Nordstrom black hole. In this case, Zerilli's 
equation can bo written as a wave equation in a complete 1 + 1 Minkowski spacetime, with a nonsingular potential 
which, being positive definite, guarantees the stability under this kind of gravitational perturbations [11]. 
When applied to the static black hole inner region r < instead, one gets a wave equation on a half 1 + 1 fiducial 
Minkowski spacetime bounded by the singularity worldline, and the potcintial has an unexpected second order pole at 
an inner point in the domain (that we call "kinematic singularity"), besides the expected divergence at the singularity. 
This makes the initial value problem for the inner region far more difficult than that for r > Tq. These technical 
difficulties, however, arc entirely analogous to those that arise in the linear perturbation problem of a negative mass 
Schwarzschild spacetime, a problem which was solved recently in [4]. As in the Schwarzschild case, the second order 
pole in the potential can be traced back to the fact that the Zerilli field $f , as defined, is a singular function of the 
perturbed metric and electromagnetic fields at the kinematic singularity (from where the name "kinematic" comes). 
Thus, an alternative field $ has to be introduced to properly analyze perturbations [4]. This is related to by an 
intertwiner operator: $ = T^f , where I = d/dx + g, and a; is a tortoise radial coordinate. In terms of $, the type 
one scalar gravitational perturbation equation is a 1 + 1 wave equation d^^/dt^ + = 0, H = d"^ /dx^ + V{x), 
with a potential V that is regular everywhere. Once an appropriate self-adjoint extension of % is chosen -and, as 
explained above, there is a unique physically motivated choice-, the evolution of initial data is unambiguously defined 
by means of an i-L mode expansion of the data. This gives a dynamics in spite of the fact that the background is non 
globally hyperbolic (see [15] for a similar approach). The intertwining technique and choice of self-adjoint extension 
is explained in detail in Section III. Previously, in Section II, we review of the basics of the Reissner-Nordstrom so- 
lution, its linear perturbations, the factorization of Zerilli's Hamiltonian, and Chandrasekhar's algebraic special modes. 

Several technical aspects of the problem are dealt with in the Appendixes. In particular, we include an algebraic 
procedure for the explicit construction of the vector and scalar zero modes considered in the paper. 

II. LINEAR PERTURBATIONS OF A REISSNER-NORDSTROM BLACK HOLE 

This section contains all the material required for the proof of instability in Section III. We first review some 
basic facts on the maximal analytic extension of the Reissner-Nordstrom solution to the Einstein-Maxwell equations 
(Section II A), and on the reduction of the linearized field equations around this solution to decoupled 1+1 wave 
equations (Section II B). Then in Section II C we calculate the perturbed Riemann tensor invariants, to determine 
the appropriate boundary conditions at the singularity for the self-consistency of the perturbation method. In Section 
II D we review from [9, 10] the factorization of the Regge- Wheeler and Zerilli hamiltonians, and its connection to 
Chandrasekhar's "algebraic special" modes, which are central in the proof of instability that follows. 
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A. The Reissner-Nordstrom spacetime and its maximal analytic extension 

The Reissner-Nordstrom space-time metric 

rfs^ = -fdt^ + ^ + ^^^^^^ ^^^^ ^ '^'^^^ d'^bdy^dy^ + r^gijdx'dx^, (1) 

is a warped product Mx-r^S'^ of a two dimensional Lorentzian "orbit" manifold times a unit two sphere. The Maxwell 
field on this space-time is 

F=9.dt^dr. (2) 
In (1), / is the norm of the Killing vector fc" = d/dt, 

/ = 1 _ ^ + ^ = ir-ro){r-n) ^ 

the latter form being useful when \Q\ < M, in which case the roots of / are positive real numbers, < < Vq, 
and correspond to the horizon radii. It is useful to keep in mind the relation between the alternative two-parameter 
descriptions of (1) 



r, = M - ^NP - Q^, ro^M+ ^ - (4) 
M = ^{n+ro), Q^ = riro. (5) 

k°' is timelike in the exterior (r > Tq) and interior < r < regions. As long as we restrict to a region where r ^ ri,ro^ 
the coordinates in (1) are appropriate. These coordinates become singular at rj and Tq, yet the Reissner-Nordstrom 
spacetime can be extended through the horizons, and new regions isometric to I: r > Tq, II: rt < r < rg and III: 
< r < Ti arise ad infinitum, giving rise to the Penrose diagram displayed in Figure 1. The stability of those regions 
isometric to III is the subject of this paper. 

Given a complete spacelike surface such as S in Figure 1, a Cauchy horizon (thicker horizon in the figure) 
develops at rj. This is the boundary of the maximal domain of development of the data, and it is connected to 
S by timelike curves of finite proper time. The solution of the Einstein-Maxwell equations is unique only up to 
the Cauchy horizon, and, although the spacetime is C°° extensible beyond it -as shown in Figure 1- the extension 
is not determined by the data in 5, and is non unique. This lack of predictability in a classical theory of fields 
moved Penrose [16] to postulate what is known as the strong cosmic censorship conjecture, according to which, for 
generic initial data in an appropriate class, the maximal domain of development is inextensible (thus guaranteeing 
the preservation of predictability). The idea that under more realistic assumptions than perfect spherical symmetry 
a Cauchy horizon would not develop, is supported by the finding that certain natural derivatives of a perturbation 
field diverge as the Cauchy horizon is approached from region II [18], and by Israel and Poisson "mass inflation" 
model [19]. An alternative, non perturbative approach was carried out by Dafermos [20]. In [20], spherical symmetric 
solutions to the Einstein-Maxwell-scalar field system are studied. The (uncharged) scalar field was added to get 
around the uniqueness Birchoff theorems in the spherically symmetric case. A characteristic problem is solved 
combining Reissner-Nordstrom data at the event horizon with generic matching data at the other null edge coming 
out the bifurcation sphere at Vg- It is shown that, generically, the Hawking mass diverges at the Cauchy horizon, and 
thus the spacetime fails to be extendible beyond it. 

The results in this paper contribute to the idea that the extended spacetime depicted in Figure 1 is an irrelevant 
solution of the Einstein-Maxwell equations. We show that a linear perturbation of the metric and Maxwell fields in 
region III, compactly supported away from the Cauchy horizon and the singularity, will grow exponentially in time, 
showing that region III is in fact an unstable static solution of the Einstein-Maxwell equations. 



B. Linearized gravity around the Reissner-Nordstrom solution 



The linearized Einstein-Maxwell equations around (2) have been analyzed by many authors, starting with the papers 
by Regge, Wheeler and Zerilli [12], generalized to higher dimensional charged black holes with constant curvature 
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FIG. 1: Penrose diagram for the ma:ximal analytic extension of the Reissner-Nordstrom space-time. All regions labeled I axe 
isometric, and so are those labeled II, and III. Regions I and III are stationary, and region I is known to be linearly stable. 
The function r is globally defined, with ro < r in I , n < r < To in II , and < r < n in III. Also r — ^ oo at the conforms! 
boundaries J^. The r = u horizon drawn thicker is a Cauchy horizon for the initial data surface S, the C°° extension beyond 
it (in particular, the two copies of region III just above it being non unique, unless analyticity of the metric is required. 



horizons by Kodama and Ishibashi [13, 14]. For polar (scalar in [13, 14], here denoted (+) following [9, 10]) and axial 
(vector in [13, 14], here denoted (— ) following [9, 10]) modes with harmonic number £ {£ = 2,3,...), the metric and 
electromagnetic perturbations of a Reissner-Nordstrom space-time can be encoded in two functions, (t, r),a = 1,2, 
that satisiy wave equations [9, 14] 



= — ^ - + =• — ^ + (6) 



with potentials 



y± = ±/3„^ + /3A' + '^/a (7) 



where k = (£ - 1)£(£ + 1)(£ + 2), /3„ = 3M + (-1)" v^]VP + 4Q2(£3]J(7Tj:2), 

= rfio. + {£-m + 2)r^ 

with / given in (3) (equation (8) corrects a typo in [2]) and x is a "tortoise" coordinate, defined by dx/dr = 1/f. Eq. 
(6) admits separation of variables $^(t,r) = exp{—iujt)^l;^{r), leading to the Schrodinger like equation 

^'i^a = nt^t (9) 

Unstable modes correspond to purely imaginary u, and thus to negative eigenvalues of the "hamiltonian" %. 



We are interested in the case M > \Q\ and < r < rj, then 

x = r+^lnf!iL^V^lnf^') 

To-n \ To J Ti-To \ n J 



(10) 
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where the integration constant was chosen so that x ranges from zero to infinity as r goes from zero to r, 
hmits 



X ~ 



3r. 



In these 



(11) 



Note form (8) that /i has a singularity at 



^9M2 + 4Q2(£- l)(^ + 2)-3M 
{l-l){i + 2) ' 



(12) 



then, from (7), we expect a singiilarity at r,; in . However, the divergences from the different terms cancel out and 
the vector potential Vi' is smooth at Vc- This is not the case for the scalar mode , which has a quadratic pole 
at rc, with a positive coefficient. The singularities at r = and r = Vc have vc;ry different origins. The first one is 
due to the spacetime singularity at this point, whereas the second one can be traced back to the definition of $f , 
which happens to be a singular function of the metric and Maxwell field first order variations at this point (see, e.g., 
(63)-(57)), this being the reason why we refer to it as a "kinematic" singularity. We should stress here that the way $^ 
is defined in terms of the perturbed metric and Maxwell fields is crucial to disentangle the linearized Einstein-Maxwell 
equations, and that Vc happens to fall outside the domain of interest r > Vg when the stability of the exterior region 
is studied. Note that Vc is a decreasing function of t and an increasing function of Q. Thus, for large enough t we 
have < Tc < rj, and V-^ is singular in the inner region, the one that we study in this paper. The "safest" mode 
is ^ = 2, for which rc > Vi, and therefore the potential regular for < r < r^, as long as Q/M < \/7/4 ~ 0.66. For 
larger Q/M rate, Vc < Vi for every harmonic mode. 

Fig. 2 depicts the i = 2 scalar potentials for some particular values of the parameters. The behaviour of the 
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FIG. 2: Scalar potentials (left) and V2' (right) for £ = 2, rj = 1, ro = 2 plotted against x. The x range was chosen in each 
case to exhibit the relevant details, beyond these ranges the behavior is that captured in equations (13) 



potentials near the spacetime singularity and the inner horizon is 



where 



The local solutions of the equation 



2 



a; ~ 



exp [- J^"^ x^ 00 



Pc. + {e-l){£ + 2)ri 



2 „/,+ 



near the inner horizon and the singularity are (for both q = 1 and a = 2) of the form 



'Acos((?)[.t1/3 J2^^q ai'^x"/3] + A sin(0)[a;2/3 J2n>o ^a;"/^] for x ~ 
&i[ exp(—A;x) + ...]+ fo2[ exp(fca;) + ...] for a; — >■ 00 



(13) 



(14) 



(15) 



(16) 



6 



where we have set a, 



(1) 





1 



-,(2) 



The differential equation (15), rewritten using r as the independent variable, has 
a regular singular point at r = 0, whose indicial equation has roots 1 and 2. The terms between square brackets 
above are just the Frobenius series solution for this equation, written in terms of x by inverting (11) (thus the powers 
of .T^/"^). The two arbitrary constants in front of them where parameterized with ^4 > and 9 G [0, 27r) for later 
convenience. Similar expansions are made for local solutions of differential equations near the singularity at different 
points below. 

Fig. 3 depicts the £ = 2 vector potentials V~ for the same parameter values as those in Fig. 2. 




FIG. 3: Vector potentials (left) and (right) for £ = 2, = l,ro = 2 plotted against x. 
The behaviour of the potentials near the spacetime singularity and the inner horizon is 



where 



V- 



c- = 



9x' 



+ ... 



C„ exp - 



('■o-r,).r 



a; ~ 



a; ^ 00 



Krj2 + l3a{ro - Ti) 



/3a + (€-l)(£ + 2)ri 



The local solutions of the equation 

H- = -k' 

which correspond to a mode oj = ±ik, are, for both a = 1, 2, of the form 

/^cos(^?)[,;:-i/3 E„>o al'^a;"/'] + ^ sin(e)[a;4/3 al^^x"/^] for a; ~ 



6i[exp(-fca;) + ...] + 62[exp(A;a;) + ...] 



for a; — )• oo 



(17) 



(18) 



(19) 



(20) 



where we have set Og^^ = 1 = Og^^ . 



C. Consistency of the linearized analysis 



In the linearized approach, a solution gab, Aa of the Einstein-Maxwell equations is replaced with a "perturbed" 
metric and electromagnetic field potential gab + shab,Aa + eBa, and the field equations are then required to hold 
to first order in e. Given that the background solution we are interested in, region III of the Reissner-Nordstrom 
spacetime, has a curvature singularity as r — >■ 0+, the perturbation treatment will certainly be inconsistent if we find 
that the first order correction to a perturbed divergent curvature scalar diverges faster than the unperturbed piece in 
the r — >■ 0"*" limit, since in this case the notion of a "uniformly small metric perturbation" is lost. This is why the first 
order piece of the Kretschmann or some other invariant is usually computed. Here we take a systematic approach 
to make sure that none of the algebraic invariant made out of the Riemman tensor acquires a correction diverging 
faster than the corresponding background metric invariant. 
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Any algebraic polynomial invariant of the Riemann tensor Rabcd can be written as a polynomial on a set of basic 
invariants, the basic invariants being generically subject to syzygies (constraints). The basic invariants are more 
compactly written in terms of the Ricci tensor Rat ■= R'^acb, the Ricci scalar R = R°'a, the trace free Ricci tensor 
Sab = Rab - QabR/^, the Weyl tensor 

2 2 

Cabcd '■= Rabcd 7;{9a\cSd]b — 9b\cSd]a) 7 TTR9a\c9d]b 

n — 2 ^ ' ^ ' n(n — 1) ^ ' 

and its dual C*^^^ := ^CabefC^'^ cd, or just using the Ricci and Weyl spinors ^j^bAb^ '^abcd- In the case of spacetimes 
with a perfect fluid or a Maxwell fleld, the basic invariants are those given in Table I (from [17].) 



Table I: Basic Riemann tensor invariants for perfect fluid or Maxwell field spacetimes 



R 




T>a 

-K 


Ri 




^ABAB^^""^^ = jS\S\ 


R2 




^^B^B^'^C^C^'^A^A = -lS\S\S-a 


Rs 




^^B^B^''c^c^''D^n^''A^A = TES\S\S'^dS''a 


Wi 




^ABCD^'''''"' = U<='abcd + iCl.JC^"^" 


W2 






mi 




*ABC/.$^^^^$^^AS = IS'^'S'^^iCacdb + iCl,J 


m2 






ms 




^abcd^Abcd^^''^^^''''^^ 


1714 




^abcd^Abcd'^^''^^'^'^'^^^^''e^e 


1715 




^ABCD^^''''^^^^^^^^'' Ab^EFEF 



In the Maxwell case, R — R2 ~ 0, and the syzygies among the remaining invariants are [17] 



R^ 



4i?3, 1714 = 0, mim2 = Rim^, m2m2m^ = Rim^mz 



(21) 



For Reissner-Nordstrom i?i and m2 do not vanish, then the syzygies imply that, as long as -Ri, Wi, W2, ^i^i and m2 
behave properly (correction does not diverge faster than unperturbed term), the same will happen to any curvature 
invariant, of any degree. Using references [12, 13], we have reconstructed the perturbed metric for each mode, and 
then calculated the perturbed invariants with the help of the grtensor symbolic manipulation package [21]. For the 
vector (axial) modes we could reduce all expressions, by repeatedly applying (6), to 



R 



Wi 
W2 

mi 

m2 



6{Q^-Mr)^ , . 6(g^ - Mr)Z ^ 



»e y«m(fc', < 



6(Q^ - Mrf 

2(Q^ - Mr)Q4 Q^Z 

^12 + 71^ ^^^C^' 

4(Q2_Mr)2Q4 A{Q'^-Mr)Q^Z 



(22) 



,,16 



le 



,.16 



where 



Z :-- 



K(/32r - 4Q2)^- 2Q{p2 + {i + 2){e-l)r)£{£+l)i)^ 



2(^2 -/3i) 



+ 



P2-P1 



(23) 
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From (20), (22) and (23), it is clear that an unstable vector mode gives an inconsistent pertmbation unless 
O'o^^ = 62 = in (20). It is only in this case that the perturbation can be uniformly bound in the whole of region III. 

For the scalar modes, the invariants can be expressed entirely in terms of the Zerilli field and its first r— derivative 
but we were not able to reduce the resulting expressions to a reasonably compact form, except for Ri, for which we 
found that 

(/^ - ) -"(«.^) - - (/^ - /¥S) -(^.«] . <-) 

where the xt is one of Chandrasekhar's algebraic special modes, introduced in the following Section (see equation 
(26)). The above formula will turn out to be very useful in the following sections. 



D. Factorization of Zerilli's Hamiltonian and algebraic special modes 



Chandrasekhar's algebraic special modes (ASM) are solutions of equations (15), (19) with real k, i.e., unstable modes 
of the perturbation equations. These modes do not satisfy appropriate boundary conditions as linear perturbations of 
region I of the Reissner-Nordstrom black hole, in agreement with the fact that this region is linearly stable. However, 
some ASM were shown in [2, 7] to satisfy appropriate boundary conditions as perturbations of the Reissner-Nordstrdm 
naked singularity, showing that this spacetime is unstable. In this Section we consider algebraic special modes as 
perturbations of region III of a Rcissncr-Nordstrom black hole, and analyze their behaviour near the singularity and 
the inner horizon. Following [9], we introduce xt defined by 



(25) 



the general solution of this equation being an irrelevant constant times 

Mrexp (^j^^ 



In terms of the functions , equation (9) reads 



Mexp 



-KX 



2^0 



(/3« + (f-l)(£ + 2)r) 



1 d^ip, 



'2j,± 



+ 



4(/3„)2 



1 d\t ^ 



(26) 



(27) 



which can easily be integrated if 



> 



(28) 



Equation (28) defines the algebraic special modes, (ASM) which give unstable solutions <I>J(t, r) = exp{kat)C,^ [r) of 
the perturbation equations, with 1^^ a solution of (27) when the term between brackets vanishes: 



Mdx 



(29) 



Note that a change of choice of R* amounts to adding a constant times x to r, and thus the constants A and B in 
(29) are unambiguously defined only after R* has been chosen. 

Alternatively, ASM can be obtained from the factorization property of the Hamiltonians 'H^ 



n+ - 

n, 



- 



(30) 
(31) 
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where 



d 



K 



A. := ^+^/^«/« + ^ ) (32) 

:= -^+(Pcfa + W^] . (33) 



dx ' V 2/3, 

^+ span de kernel of AaBa, with x+ in ker Ba, whereas (^^ span de kernel of BaAa, with x~ in ker ^q,. 



1. Algebraic special vector modes 



The asymptotic behaviour of Xa ^^^^ the spacetime singularity and the inner horizon is 



X 



-1/3 



Xo 



(3r,ro)i/3 

f + -!)(£ + 2)ri]exp(-^) 



The vector r modes are 



/3c + (^-l)(^ + 2)r 
^exp(^) 



/■X r^exp^- 

Jr, + - 1 



a; ~ 

X — >■ 00 

da; 



(34) 



' + 2)r 



2 ■ 



(35) 



Since /3i < and < /32, for a = 2 the integral in (35) near x = oo (r = rj) diverges, we can give R* any /iniie value, 
and the asymptotic behaviour will depend on whether i?* = or not: 



'a=2,fl»=0 



'^4/3 

exp 



( KX \ 
2/32 j 



a; ~ 
a; ^ 00 



(36) 



'a:-i/3 
exp 



/ KX \ 



X 
X 





> oo 



(37) 



Since those type-2 vector ASM that grow slower than as r 0+ blow up at the inner horizon, it follows from 
the analysis of invariants in the previous subsection (equations (22)- (23)), that a pure mode $2 = 6*^^* ^2" = 
e*^^* [A2 X2 +-^2^ '''2R'] cannot be consistently treated as a first order perturbation on the entire Reissner-Nordstrom 
region III. 

If a = 1 the integral in (35) near infinity converges, thus, we can give R* any finite value, or take R* = 00. The 
asymptotic behaviour will be 



r4/3 



exp 



( —KX \ 
2h) 



+ ■ 



a; ~ 
a; ^ 00 



(38) 



'a=l,ii*#0 



-1/3 



exp 



(_ ^x 1 
20 J 



+ 



a; ~ 
a; ^ 00 



(39) 



-1/3 



exp 



( KX \ 
201 ) 



+ 



a; ~ 
a; 00 



(40) 



It follows again from equations (22)-(23) that a pure AS mode $j = e'^^* Ci = e'^'* [^1 Xi + -^i r*] cannot be 
consistently treated as a first order perturbation in the entire Reissner-Nordstrom region III. 
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Note that some care is required when interpreting equation (35) for the a = 1 vector mode, due to the singularity 
of the integrand at Tg. Take, e.g., the case R* = oo. The one form under the integral sign in 



'l,oo 



{rc - r) 



(^-l)(£ + 2)r exp(^) 



7 — ~ — ^2 *^^P ( ^ ) 
{r-Tcy \/3i 



(41) 



can be written as 



+ + Z{r) dr, Z{r) the regular function obtained by subtracting the second and first 
order poles. The integration constants at both sides of rc can then be adjusted such that 



A-B 



''"1,00 



1 1 ir- is well defined across Tc 



(r - rc) In \r - rd + (^B In |r,; - - 7-37;) {r - rc) + (r - rc) J^' Z{r)dr 



(^-l)(£ + 2)r exp(^) 



2. Algebraic special scalar modes 
The asymptotic behaviour of xt near the spacetime singularity and the inner horizon is 



M 



a — 



a; ~ 
X ^ 00 



The asymptotic behaviour of r^jj. depends on the choice of R* in 



r exp 



I KX \ 



Po. + {£-l){£ + 2)r 



" [Po, + {e-i){e + 2)r]^dx 



exp 



( KX \ 
0j 



If a = 1 the integral (43) diverges near infinity, then R* is restricted to finite values, and 

fxVs + ... 

exp I 



'a=l,fl»=0 



(KX \ 



a; ~ 
a; — >■ 00 



(42) 



(43) 



(44) 



If a = 2 there are three possibilities: 



'a=2,fl»=0 



r.1/3 



exp 



201 J 



+ ■ 



I exp 



( KX \ 
2/32 j 



a; ~ 
a; — >■ 00 



X ~ 
a; — )■ 00 



(45) 



(46) 



'^1/3 

exp 



( KX \ 
2/32 J 



a; ~ 
a; ^ 00 



(47) 



'^a=2,iJ*=oo 



a/3 



exp 



(_ KX \ 
2/32 J 



+ 



a; ~ 
a; — )• 00 



(48) 



The only scalar AS modes that behave appropriately near the inner horizon are xti=i ^^^^ "^0=2 b*=oo- Since the 
integrals defining the latter one are non elementary, we proceed with xj=i) fo'^ which an explicit reconstruction of the 
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perturbed metric and invariants is relatively simple. No invariant that is trivial at order zeroth develops a first order 
correction. The non zero invariants to first order for the xt perturbation are 



mi 
m2 



6(Mr - g2)2 



2\3 



Q{Mr-Q ^) 



+ e 



3M^2^(^+1) {Mr-Q^) 

2/3ir7 

9M/32 {Mr - g2) 



11 



APir 
MQ%^(^+1) 



4/3ir 



11 



4(Mr - Q2)2q4 
ms = ^ e 



r.16 



(t-x) 



/3ir 
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(49) 
(50) 

(51) 

(52) 

(53) 



This algebraic special mode satisfies all our requirements for a self consistent first order formalism in this singular 
spacetime. Thus, we will restrict to type one scalar perturbations from now on. Xi is an example of a spatially 
uniformly bounded perturbation that grows exponentially in time, and thus a signal of a gravitational instability. In 
the following Section, we will show that this mode can be excited by a generic perturbation that is initially compactly 
supported within region III. The treatment will follow closely the case of the Schwarzschild naked singularity treated 
in [4]. 



III. PROOF OF THE LINEAR INSTABILITY OF THE INNER STATIC REGION 



The linear instability of a static spacetime is established once an unstable mode is found. We have shown in the 
previous section that, out of the four modes existing for every harmonic pair (£, m), the type=l scalar mode admits an 
unstable solution to the linearized Einstein-Maxwell equations -Chandrasekhar's AS mode- that can consistently be 
treated to first order in the whole domain of region III, < r < r^. The purpose of this section is to show how generic 
initial data with compact support in region III excites this mode. Although this problem is trivial for perturbations 
in region I, it exhibits a number of unexpected difficulties when dealing with perturbations of region III. Zerilli's wave 
equation for this mode: 



= — + y, $1 



0*$ + 



(54) 



has a potential with a singularity at Tc given in equation (12), that generically falls in region III (see left panel in 
Figure 2). The origin of this singularity (a second order pole) in can be traced back to the definition of $f 
[12, 13]. The first order variation of the electromagnetic field has 



_ dAdYe^ , 



dr de 



then A must be smooth for 5F to be smooth. Since 



fe(r-re) + 
8Qr ^1^*'"^^' 



(55) 
(56) 

(57) 



we conclude that, for generic smooth perturbations, has a first order pole at r = Tc, and can be Laurent expanded 
around rr as 



E 

k>-l 



Ck{r-rcf 



(58) 



The variation of gtr can be simplified to the form [12, 13] 



Y. ^ 
at 



r'-^+B^t 
or 



(59) 
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where 

{{i - 1) + 2) - 3 M (-3 + e^ + e)r^ + (2 {£ + 2) _ 12 M^) + 13 Q'^Mr - 4 Q*) 



(r2 -2Mr + Q2)((^_ 1) (^ + 2) +6Mr-4g2) 



rV9M2 + 4Q2(£_l)(^ + 2) 
((^-l)(^ + 2)r2+6Mr-4Q2) ^"^^^ 

It can be checked that the (r — rc)~^ coefficient of the series expansion of (59) vanishes. The (r — rc)~^ coefficient 
will vanish if 

CO _ (^-l)^(^ + 2)^(2M-^i) 

c_i pi{2{l-l){£ + 2)M+{P+l + 2)l3i) ^ ' 

It turns out that (58) and (61) are not only necessary but also sufficient conditions for the perturbations of the metric 
and electromagnetic fields to be smooth. Thus, we have proved that the Zerilli functions <I>^ corresponding to smooth 
type-1 scalar perturbations are those admitting, at any fixed time, a Laurent expansion (58) around Vc satisfying the 
condition (61). 

We also need to check that the perturbation does not change the character of the singularity. As explained in the 
previous section, this guarantees the self-consistency of the linearized theory. It then follows from equation (24) that 
we must demand that, for some positive 5 and N 



dr dr 



<N, if 0<r<5 (62) 



We will show in the following Section that this condition is also sufficient to assure that all the invariants behave 
properly. Thus, we arrive at: 



Lemma 1 In order that the metric and electromagnetic scalar type-1 perturbations be smooth and the linearized 
approach be self consistent, the Zerilli function has to satisfy (58) and (61) whenever Vc < ri. It also has to satisfy 
condition (62), and decay properly as r ^ r^. In particular, both initial data functions ^i{t = 0, r),$^(t = 0,r) 
must satisfy all these conditions. 

The rather odd initial value problem that these conditions pose is in fact very similar to that of the propagation 
of scalar gravitational perturbations on a negative mass Schwarzschild background, which also has a "kinematic" 
singularity. This latter problem was worked out in [4] using an intertwiner operator [24]. We will apply the same 
technique in what follows. 



A. Basics of the intertwining technique 



Consider a wave equation with a time independent potential: 

on a domain t G R, x G {a, b) where a = —00 and/or 6 = 00 is a possibility. An intertwinner for this equation has the 
form [4, 24] 

with tpx satisfying 

n^x = Exyji. (65) 

The above equations neither assume that H is self adjoint in L^(^[a, 6), dx) nor that ipi is an eigenfunction of such an 
operator, ipx in equation (65) is just any solution of this differential equation, without any consideration on boundary 
conditions, boundedness or finiteness of some norm. 
In terms of 



(66) 
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equation (63) reads 

= ^+U^ 

^:=-^ + V^ (67) 
V := 

ax 

That is, if $ is a solution of (63) with initial data 

mt = 0,x),^t = 0,x)), (68) 
then 1$ =: $ is a solution of (67) with initial data 

($(t = 0,x),^(t = 0,x)) = (X$(t = 0,x),T$(t = 0,x)). (69) 

The general idea of the intertwiner technique is to use this fact to search for an appropriate intertwiner such that V 
is simpler than the potential in the original problem. 

The operator X has a nontrivial kernel spanned by i/jx, so information is lost when switching from $ to $ := 1$. 
The operator 

X=^+., (70) 

can easily be seen to map solutions of (67) onto solutions of (63). A straightforward computation shows that 

II = Ex-n (71) 

Since I and I have non trivial kernels, information is lost when applying these operators. However, in the case of a 
solution of the wave equation (71) implies that 

Jl> = J J$ = {Ex - H)$ =Ex^ + ^, (72) 

Thus the information lost is precisely the $ initial data. The above equation can be regarded as a ODE on $ for 
every .X, and can easily be integrated to give $ back: 

If Ex > 0, 

$(i, x) = /* sin (^/Ei{t - t')) l^t', x)dt' + cos(/^i)$(0, x) + ^'"^^^^ ^(0, a;); (73) 

V Ex Jo ^ ■' V Ex 

if Ex < 0, 

^{t,x) = -jl= l\mh(^/^{t-t')) U{t', x)dt' + cosh(v/^0^(0. x) + '^"^'^;^^^ ^$(0, x); (74) 
y—Ex Jo ^ ' y—Ex 

if Ex = 0, 

$(t, x)= j^{t- t') X^t', x)dt' + $(0, x) + 1 $(0, x) = (^J^ ±^{t", x)dt''^ dt' + 1^{0, x) + $(0, a;). (75) 

We conclude that we can solve equation (63) subject to (68) by means of the following procedure: 

1. From the initial $ data (68) construct initial $ data using (69). 

2. Find the solution $ of equation (67) with initial data (69). 

3. Apply (73)-(75) to obtain the solution $ to the original equation. 

Note that the evolution problem is solved in step 2, and that the initial $ data is used twice: in steps 1 and 3. 
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B. The initial value problem for perturbations of region III 



Why would be one be interested in solving (63) using the complicated intertwiner method? The intertwiner is 
certainly useful if V is simpler than V. Our motivation, however, comes from a deeper problem: there is no available 
theory to deal with the initial value problem possed in Lemma 1. Unless rc> U, which may only happen for a finite 
number of harmonic numbers £, the function space in Lemma 1 is unrelated to any recognizable Hilbert space, and the 
Zcrilli wave equation has singular coefficients in region 111 (on top of this, there is the issue of non global hypcrbolicity 
of region 111, even when rc > rj.) A similar situation is found when studying perturbations of a negative mass 
Schwarzschild spacetime. In this case, the problem was solved using intertwiners [4]. Less sophisticated approaches, 
such as using algebraic redefinitions of the variables, can be shown to fail. 

As we show in Appendix A, for any (even complex) Ex, the intertwiner (64)-(65) produces a V that is smooth at 
r = Tc, while sending functions satisfying (58) and (61) onto functions which arc smooth at r — r^- Also, there are 
options for il>i for which I sends initial data (as characterized in Lemma 1) onto a subspace T) C ^^((0, oo), dx) where 
"H is a self-adjoint operator. This allows us to use the resolution of the identity for T-L to solve the hat wave equation 
by separation of variables, by expanding the initial data using normalized eigenfunctions of %, 'HiIje = E'tpE- 

l>(t,x) = Y,"-Eit,x)i^E{x) (76) 

E 

mx) = Y^a%i>E{x) (77) 

E 

Ht,x) = J2^eMx)- (78) 

E 

Here the coefficients are obtained by integrating against the complex conjugate of '^e{x), and the wave equation then 
reduces to an infinite set of ODEs: 

CLE = -EaE 

aB(0) = d%:= j tpE Ht = 0,x)dx (79) 
aE{0) = a%:= J ^pE ^it = 0,x)dx. 

whose solution is 

I cos 

aE{t) = < a% + 1 d% ,E = (80) 

[a% cosh{y/^t) + d°E {-E)-^/'^ sinh(^/^^) ,E <0 

The above equations define the evolution of the fields outside the domain of dependence of the initial data. This 
same technique was applied, e.g., in [4, 15], in similar contexts. A subtle issue is that of defining the domain 

T) C L2((0, oo), dx) where "H is self-adjoint. This problem is identical to that of quantum mechanics on a half axis, 
treated in detail in the first reference in [22] . Consider the two dimensional vector space of local (Frobcnius) solutions 
of the eigenvalue equation T-Lijj = Eijj, ip ^ 0, near x = 0. Given that an overall factor on is irrelevant, the space of 
local solutions can be regarded as a circle (this is why we used A cos(^) and A sin(^) for the two arbitrary constants 
in equations (16), (20), etc. 9 E [0,27r) labels points in this circle.) 

If any eigenfunction is square integrable near x = we say, following [22], that H belongs to the "limit circle case". 
In this case, H will be a self- adjoint operator only after restricting to a subspace Vg^ C £^((0, oo), dx) of functions 
behaving near x = as local eigenfunctions with a fixed 9 = 9o- Equations (76)-(80) will then hold in Vg^, for initial 
data in this space. Note, however, that initial data of compact support belongs to Vg for any 9, and evolve in a 
different way if aomc ^ 9o is chosen in (76)-(80). Of course, the solution will be different only outside the domain 
of dependence of the initial data, but still there is an ambiguity, which must be resolved. Physical input must then 
dictate what the right choice of 9 is in order to get rid of this ambiguity. 

The other possibility is that the hamiltonian piece of the wave equation belongs to the "limit point case", i.e., that 
there is a single 9 value giving local solutions which are square integrable near x = 0. In this case we say that H is 
"essentially self-adjoint" (since it is only self-adjoint in the domain defined by this particular 9 value) and there is no 
ambiguity in the dynamics. This would be the case if one of the roots of the indicial equation of of the Probenius 
local solution of the hamiltonian eigenfunction were less than —1/2. 
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For the scalar gravitational perturbation problem, the situation is that of a limit circle hamiltonian. The self 
consistency condition (62), however, singles out a unique 9. With this choice, the degree of divergency as x — >■ 0"*" 
not only of but also of all the remaining algebraic invariants of the Riemann tensor, gets controlled, and, since 
the evolution (76)-(80) preserves the local behaviour at a; = 0, the invariants will stay properly bounded near the 
singularity at later times. The dynamics is thus unambiguous once we enforce the self consistency condition (62). 

In [4], a "zero mode" {Ex = 0) was used to construct the intertwiner and produce a self-adjoint %. The resulting 
hamiltonian has a negative energy eigenvalue, and thus exhibits the instability of the spacetime. The Zerilli field 
is then recovered using equation (74), from where it is clear that the exponential growing in time of $ shows up 
in the metric and electromagnetic field perturbations. We have tried this same approach here, and found that an 
appropriate zero mode can be explicitly constructed for ^ = 2 and, as happens in the Schwarzschild naked singularity 
case. i-L has a smooth potential and contains a negative energy eigenvalue, at least for some Q/M values for which 
I'c < fi (see Section IIIC). Since generic perturbation initial data with compact support in region III will have a 
nonzero projection onto the £ = 2 type-one scalar mode, this is certainly enough to show that such initial data will 
excite unstable mode in these cases. 

However, we were not able to prove that for arbitrary i and Q/M such that fc < r, there is a zero mode intertwiner 
which does not introduce a new singularity in V (although it is still trivial to show that any intertwiner washes out the 
singularity at r = Tc, see Appendix A.) A new singularity would be introduced if tpx had a zero for r e (0, fg) U (fg, r^). 

For this reason, in Section III D we exhibit an alternative intertwiner for which computations can be carried out 
explicitly for any i. This uses ipi = xf > Chandrasekhar's mode, and gives, for any Q/M and i, H = Hi (the 
Hamiltonian for type-1 vector perturbationsl). Since H is positive definite in this case, the hat wave equation is 
stable, and the scalar instability shows up only when reconstructing the Zerilli field using (74). This is so because 
those factors inside the integral which are exponential in t do not cancel the exponential factors outside the integral 
(as it would happen if the original wave equation were stable). These factors will then show up in the metric and 
electromagnetic field perturbations, the Riemann tensor, and its invariants. 



C. The £ = 2 zero mode intertwiner 



In [4], a "zero mode" (solution of Hip = that is not necessarily normalizable or well behaved at the boundaries) was 
used to construct an intertwiner to deal with the initial value problem for the scalar mode negative mass Schwarzschild 
perturbations, which has difficulties similar to those found in the present case. 

We may try the same approach here, however, given the complexity of the potential in Hf , it is rather difiicult to 
obtain the solutions of 



(81) 



for tjjx = ijjo required to construct the Ex = intertwiner X (64). One possibility is to use the relations (30)-(31) to 
obtain scalar zero modes from vector ones, since, for a vector zero mode, 



(82) 



it follows from (30)-(31) that Aitjj^ is a scalar zero mode. Since the vector potential V{' is much simpler than V-^ , 
there is some hope that we could carry on calculations in a more explicit way using this idea. This is indeed the case 
for £ = 2, for which the general solution of equation (82) (see Appendix B) can be shown to be 



= ^cos(a) <» + ^(Q2 _ _ 9l 



(a) |r 



Q' 



In 



-r + M 



- 1 - In 



-r + M 



+ 1 



. 9 12M3-2(^2-/3i)(M2-g2) 
12 + 3 (/3i -2M)r+ (3M/3i - 2{M^ + 2Q^)) + ^— — ^ ^ 



• (83) 



where A and a are arbitrary constants. Note that the £ = 2 intertwiner operator constructed using V'l = i>o = 

in (64) will depend on a but certainly not on A. It can be easily shown, however, that V is smooth at r = rc 

irrespective of the choice of a, the first and second order poles in V^^ being canceled by the poles in dg/dx (Figure 4). 
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This, of course, is to be expected from the more general considerations in Appendix A. The asymptotic behaviour of 
V near the singularity and the inner horizon is also independent of a: 

fA + ... ,a;~0 

V-in ( (r^-r.)A (84) 
I^Cexp (^- ^ "^a'-" j ,a;->oo 

Near the singularity the eigenfunctions of H behave as 

Vi = Acos(^)[x-i/3 a(,i)x"/3] + A sin(^)[x^/3 ^ a^f)^"/^] (85) 

n>0 n>0 

Thus a belongs to the limit circle case, and only restricting to a subspace Vg^ of functions behaving as (85) with a 
fixed — 9o value does H becomes self-adjoint. 

We will make the choice 6 = n/2 of slowest decaying functions. This condition is certainly preserved by the hat wave 
equation (see (76)), and implies 

^$ = ^afcr^ a2 = ^au (86) 
fc>i ^ 



near r = 0. If $^(0, r), ^^(0, r) also admit expansions like those in (86), then, using (74) follows that 

A 

4Q 



$+(^,r) = ^afc(^)r^ a^it) = ^ a,{t) , (87) 



k>l 

for all t, a condition that can be shown to guarantee that all algebraic invariants of the Riemann tensor behave 
properly near the singularity. 

Regarding the choice of a in (83), although the results do not depend on the intertwiner that we use, it is certainly 
easier to understand how the instability is excited by an initially compactly bounded perturbation if we use 

tan(a) = ^ ^ (33) 



3Q'/32 In (^^^3) + 2(3M^2 - 16(M2 - Q^))^M^-Q^^ 



since in this case the resulting intertwiner will send xt onto the the Hilbert space 15^/2 selected by the self consistency 
argument, and thus Ixt '^il^ one of the eigenfunctions of ThL (it will actually be the only negative energy H 
eigenfunction). The transformed potential V, together with Ixi for this choice are given in Figure 4 below for some 
specific Q and M values. An explicit expression for V can be readily obtained using equations (64), (67), ipo = -^iV'o ) 
(83) and (88). 

Now suppose some perturbation data {^i{t = 0,x),^f{t = 0,x)) of compact support is given. The hat wave 
equation data (T$^(t = 0,x),T<i>^(i = 0,x)) will be of compact support and then it will belong to D-^ji- Expanding 
it using (77)-(78) will generically give a nonzero projection onto the fundamental, unstable ?i eigenfunction Xxf> 
and thus, from (80), an exponentially growing term in (76), which survives when $f is reconstructed using (75) and 
shows up in the metric and electromagnetic field perturbations. 



The use of a zero mode intertwiner has some drawbacks: although we can show that the kinematic singularity is 
absent from Y (see Appendix B), we do not have a complete proof, even for ^ = 2, that the zero mode '01 = "0^ hs-s 
no zeroes in (0,rc) U (rc,rj), which would introduce new singularities in V . However, for i! = 2, we have numerically 
verified that this is the case for a range of values oiQ/M. A particular example of a smooth F for ^ = 2 is that given 
in Figure 4. 

In the following Section, we show that all these issues can be avoided by using an alternative intertwiner that allows 
explicit calculations for every harmonic number and charge and mass values. 
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FIG. 4: £ — 2 potential V -continuous line- for the transformed Zerilli equation (67), and the transformed of the unstable mode 

xt ■ In this example rt = 1 and To = 2. 

D. Intertwining using Chandrasekhar's algebraic special mode 

We can apply the intertwining technique using Chandrasekhar's algebraic special mode Xi in (64), for which 

\ 2 

K 



Using this intertwiner has a number of advantages. Contrary to what happens for the zero mode, we have an explicit 
expression for foi' (Hx^ry harmonic number i, equation (26). It is also simple to construct V using (7) and (25) (a 
prime denotes derivative with respect to x): 

f I Dili , i\„2 o „ , /i/o2\ 




= = ^ + \y - /3ir + 4Q2) (90) 

(This relation between the scalar and vector modes was first noticed by Chandrasekhar [9, 10].) The fact that 
V = , the type-1 vector potential, simplifies the analysis considerably, since it is clear that is smooth in (0, rj) 
for any value oiQ/M. 

From the first line in (20) follows that i-L = belongs to the limit circle case. As explained above, a choice Oo has 
to be made to fix the domain Vg^ C L^((0, oo), dx) where H is self adjoint. However, for type-1 scalar perturbation 
the consistency requirement (62) (see also Lemma 1) reduces to: 

l^t <N, if < r < 5 (91) 

which rules the x^^^"^ piece of (20), and forces 9 = tt/2. Once again, this condition is preserved by the hat wave 
equation (see (76)), and imply (86) and (87) for initial data satisfying this condition (where now X has to be 
understood as the intertwiner made using Chandrasekhar's mode) Condition (87) guarantees that all algebraic 
invariants of the Riemann tensor behave properly near the singularity at later times. 

Since V (and thus H.) is positive definite, the hat wave equation is stable, and its modes oscillate in time. The insta- 
bility of arises as a consequence of the fact that, generically, the exponential terms in the integrand of (74) do not 
cancel out with those outside the integral. As a trivial example, take ^i{t = 0, x) and ^i{t = 0, x) both proportional 
to Xi (this certainly passes all the requirements in Lemma 1). Since Ixt = 0, the initial data in hat space is trivial, 
then $(f, a;) = for all t, and (74) reduces to the last two terms, which are generically exponentially growing for large t. 

A final observation is the nontrivial fact that, unlike the Zerilli field, the intertwined variable has a geometrical 
significance, as, from (24), it gives the first order variation of the Riemann invariant -Ri: 

^^t = ^ = -^SR, (92) 



18 



This gives further significance to the dual relation between vector and scalar modes first found by Chandrasekhar, 
which was limited to some observations on their mode spectra. Equations (24) and (90) prove that the field giving 
the first order variation of Ri (times Q'^/r^) associated to a scalar mode perturbation is a solution of the vector mode 
perturbation equation of the same harmonic number. 

IV. CONCLUSIONS 

We have proved that the inner static region < r < of a Reissner-Nordstrom black hole is unstable under 
linear perturbations of the metric and electromagnetic field. More precisely, we have shown that a perturbation with 
compact support within this region will excite unstable type-1 polar modes, of which there is one for every harmonic 
number (i?, m). This instability is relevant to the strong cosmic censorship conjecture, according to which this region of 
the black hole, lying beyond the Cauchy horizon (of a Cauchy surface like the one in Figure 1), should be disregarded, 
as it could not arise as a result of the collapse of ordinary matter departing from spherical symmetry. 
This result has implications on some simple models of halted collapse of a pressure- less charged perfect fluid star [25] , 
according to which the world tube of the surface of the star traces a path going from the (right copy of) region I in 
Figure 1, through region II, into the left copy of region III, upper copy of region II, then upper right copy of region 
I. To the right of this curve, this spacetime agrees with the extended Reissner-Nordstrom spacetime, which contains 
an entire copy of region III, thus being unstable. We are currently studying these models in more detail. 
The difficulty in establishing in a rigorous way the Reissner-Nordstrom instability lies in the fact that the field variable 
which succeeds in disentangling the linearized equations, the Zerilli field $f , happens to be a singular function of 
the perturbed fields in region III. This cannot be cured by any simple field redefinition, but requires the use of an 
intertwiner operator I = d / dx + g that maps to a smooth field $ := 1$^. The information lost due to the nontrivial 
kernel of X is entirely contained in the initial data, and thus is available. The evolution of perturbations on the non 
globally hyperbolic background is well defined by using the spectral theorem and a unique self-adjoint extension of the 
spatial piece of the wave operator that gives the dynamics of the $ field. We should comment here that intertwiners 
in the context of linear perturbations were first considered in [24] , while they were first used to deal with the issue of 
the singularities of the Zerilli field in the proof of the instability of the Schwarzschild naked singularity in [4]. The 
idea of defining dynamics on non globally hyperbolic backgrounds by using a suitable self adjoint extension of the 
spatial piece of the wave equation together with the spectral theorem was first suggested in [15]. The main difference 
between the cases considcircid in [15], and the Reissner-Nordstrom and negative mass Schwarzschild cases, lies in the 
fact that the spatial operators ( "Hamiltonian" ) in the last two cases are not positive definite. A difference between 
the negative mass Schwarzschild and the Reissner-Nordstrom cases, is that the intertwiner used in the first case gives 
a Hamiltonian with a single self-adjoint extension (limit point case in [22]), whereas the one for Reissner-Nordstrom 
corresponds to the limit circle case in [22]. 

Two different intertwiners were used: one constructed out of a zero mode, the other using one of Chandrasekhar's 

algebraic special modes. The first intertwiner has the advantage of exhibiting the instability in a rather obvious 
way, and the drawback that we lack explicit expressions for the interwtined potential, or a proof of its smoothness 
within the relevant parameter range. The intertwiner that uses Chandrasekhar's "hides" the instability, which is made 
explicit in the metric reconstruction process through the original Zerilli field. This mode allows explicit calculations 
for every harmonic number and Q and M values. It also exhibits a very interesting connection between vector and 
scalar modes: the intertwined field 3> gives the first order variation of the Riemann invariant Ri (see equation (92). 
An alternative way of stating this is that the first order variation 5R\ of R\ under scalar perturbations is a solution 
of the Zerilli vector perturbation equation. 

The results presented here adapt easily to the case of a super-extreme (|Q| > M) Reissner-Nordstrom spacetime, and 
thus can be used to fill in the details left untreated in [2] to show that a perturbation of an overcharged Reissner- 
Nordstrom spacetime, compactly supported away from the singularity, will excite modes that grow exponentially in 
time. This, of course, is relevant to weak cosmic censorship, this being the original motivation for our work. 
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Appendix A: Intertwined potential 



In this section we analyze the local behaviour of the intertwined potential and 'ip given in equations (64), (65) and 
(67) as the singularity, inner horizon, and r = Tc are approached. This is done by studying the behaviour of generic 
solutions of the equation 



+ V+tPx = Ei^t 



for r ~ 0, r ~ Tc, and r ~ri. 

Note that may be written as, 

/ 



showing explicitly the double pole at. 



dr 



r = rr, = - 



(A + (^ - m + 2)r) - 2/(r)/3i(^ - 1)(^ + 2) 



/3i 



(Al) 



(A2) 



(A3) 



1. Behaviour of V 

Near r = Tc, the general solution of (Al) is of the form, 

(^-l)2(^ + 2)2(2M-/3i)co 



1 



r -rc 



Co 



I3i {2{e -!){£ + 2)M +{e^+£ + 2)^i) 
SEx^j CO 



{2{£ - + 2)M + (£2 £ + 2)/3i)- 



r(r" - Vcf + C3(r - rc)^ + C4(r - Tc)^ + . . . 



(A4) 



where cq and C3 are arbitrary constants, and C4 and higher coefficients in the series are linear combinations of cq 
and C3 with coefficients that depend on Ex, Q, M and L Note that the generic local eigenfunction above satisfies the 
requirement (61). 



If we use the generic V'l given above to construct the potential V, 



V = V~2f- 



d 



I d 



^2 



dr \ ipx dr 

a straightforward computation shows that, provided qq ^0, near r = Tc, 

8fc2^f + + 2f{l - + £ + 4)/32 - 20M/3i - 12(^ + 2){l- 1)M2] 



V = 



+ 0{r-re). 



(A5) 



(A6) 



This means that, for any Ex, an arbitrary solution of (Al) with ao 7^ gives an intertwined potential V that is 
smooth at r = rc- We notice also that, as can be checked, if ao = 0, the second term in the R.H.S. of (A5) does not 
compensate the double pole in V, so that V is also singular. 

We consider next the local behaviour near r = 0. The general solution of (Al) admits an expansion in powers of r 
of the form. 



V'i(r) = air + a2r^ + 



(f-l)(£ + 2)((£-l)(£ + 2)Q2 + M/3i) M 



+ aAr^ 



(A7) 



where ai , and a2 arc arbitrary constants and the higher order coefficients depend linearly on them. A dependence on 
Ex appears first at order r^ . This result implies that, near r = 0, assuming a\ ^ 0, we have. 



- _ 4 _2 2 / M 2a: 

9^ " 3V3 ai J 



-5/3 



(A8) 
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while, if ai = 0, 



10 



32/3 / 10e(£ + 1) - 12) M(28M - 5/3i) 



)2/3 



+ ■ 



j;-4/3 + 0(^-1) 



(A9) 



Finally we consider the behaviour near r = r^. The cases Ex ^ and Ex = require separate treatment. For 
Ex ^ 0, since the potential vanishes for r = rt, the leading terms of the two linearly independent parts of the solution 
for real Ex are of the form, 



ipx = Ci {Vi - r) 



-rlVEi 



= Ci exp (^^/E^x^ + C2 exp (^-y/s^x^ 



(AlO) 



where Ci, C2, Ci, and C2 are constants. For Ex = 0, on the other hand, the general solution admits an expansion of 
the form, 



tptir) = ao + 



{{n - ro)Pi + ktD ao - n (/3i 



2)(£ - 1)((2£2 + 2£ - l)ri + 2ro)) 60 



+a'2{r - nf + . . . 



+ ln(ri — r) 



60 



r-i(n-ro)(/3i+n(^+2)(£-l)) 
((r-i - ro)/3i + Krf) 60 



(All) 



ri(ri-r„)(A+n(^ + 2)(^-l)) 



(r - Tj) + 62(r - Tj) +... 



where 02, &2 and higher order coefficients depend linearly on ao and 60 • 
The Ex = transformed potential behaves as 

where the dots indicate terms that vanish as a (r^ — r). In terms of x, for bo ^ 0, this implies, 

V{x) = ^ + ... 



(A12) 



(A13) 



2. Behaviour of ij) 

We consider now the behaviour of the intertwined field -tjj (64), which can be written as 

^ = T^ = f^4{^). (A14) 
We will use the following result, whose proof is straightforward: 

Lemma 2: Assume that tp and ipx admit a Laurent expansion 

^ = (r - roY 5^ afe (r - r„)^ = (r - Y.'^kir- Vo)'' (A15) 

fc>0 fc>0 

where ao = 1 = , and p is any integer number. If s is the highest number for which Uk = for every k < s {s 
measures the degree of contact of these functions at Tq, and, generically, s = 0), then 

= (r - ro f Yldkir - r„)^ 4 = s(aa+i - aj+i). (A16) 

Consider first the action of the intertwiner on a function ip satisfying (58) and (61). Since, as follows from (A4), 
■07 satisfies this same condition, generically tjj and ijjx have (as functions of r) degree of contact s = 2, in the notation 

of Lemma 2, and thus ip = /V'x^ (^) ^^^^ smooth at r = Vc- 

The local solutions of Hiip = i?V-'j are of the form ■0 = X)fe>i '^fe^'^ with 02 and 04 arbitrary, Ofc independent of E up 
to = 7. Consider the action of I on a function like this further subject to the condition (86). If the intertwiner also 
satisfies (86), the degree of contact will be s = 5, then ■01^ (^) will be ©(r^), and thus 'tp will be ©(r'') = 0(x*/^). 
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Appendix B: The scalar and vector zero modes. 



In this Appendix we describe a procedure that allows the construction of the zero mode solutions for both vector 
and scalar modes. We start with the vector zero modes by first considering the differential equation they satisfy. In 
accordance with (6, 7), and (9) with w = 0, this is given by, 

We notice that (Bl) has regular singular points for r = 0, r = rj and r = To, and no other singularity. From now 
on we will write simply -0 for ipQ . 

A simple analysis of the indicial equation shows that near r = this equation has two independent solutions, one 
behaving as r~^, and the other as r^, both admitting a power series expansion. We consider therefore an expansion 
for ipQ (r) of the form, 



oo 



Replacing in (Bl) we find that we must set 02 = 0, and, 

(32 £{i+ 1)132 £(£2 _!)(£ + 2) 
= "4Q^^° ' = 24Q^ ' = 24Qi "° ^^^^ 

The coefficient 05 can be chosen arbitrarily, in accordance with the indicial equation, and is given by, 

^2 - i ){e2 -4){£ + 3) _ /3i - 1 6M 
T44Q6 

For the remaining coefficients, including oe, we find a three term recursion relation of the form, 

-{£ + J - 2){i - J + 3)a,_i + (A - 2(j - l)(j - 3)M)a, + (j + l)(j - A)Q^a^+, = (B5) 



«6 7777^6 S752 1^4) 



and, therefore, all the coefficients are determined once ao, and 05 are given. But this implies that for any given and 
ao ^0, we may choose 05 such that = 0, and then all coefficients for j > £ + 3 vanish. Calling V'a this solution 
we have, 

Mr) = -Veir) (B6) 
r 

where 'Pe{r) is a polynomial of order £ + 2. The lowest order polynomials are. 



r5 (B7) 



P2 ^ 

r + 


P2 3 


1 


4Q4^ 




P2 , 


/32 3 


5 




2Q4'' 






5^2 3 _ 


15 


6Q^ 





{p2 + 10M)Q 

r,(r) = 1 - M — "r^ I 21(24M + /32) 5 ^2 



(B8) 



where we have fixed, for simplicity, 'P^(r = 0) = 1. With this normalization, the polynomials are positive and 
decreasing functions of r near r = 0. We shall now prove that they have no zeros in the interval < r < rj. We write 
(Bl) in the form, 

^ {2nro - (r, + ro)r) dj^g ^ {£{£ + l)r^ - /3ir + AnVp) ^_ 
dr'^ r{ro — r){ri — r) dr r'^{ro — r){ri — r) " 

and notice that ■^a can have only simple zeros in < r < because the coefficients in (B9) are regular functions of r 
in < r < Tj. Since sufficiently near r = we have ^pa > ^ and dtpa/dr < 0, and the coefficients on the RHS in (B9) 

are both positive, the sign of d'^ipa/dr'^ is not fixed. We notice however that at the first zero of V'a for r > we must 
have dtpa/dr < 0, and therefore, we also have d'^tpa/dr^ < 0. Since to the right of such a zero, and as long as r < r^. 
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(since the coefncients arc still positive) we must have both dtpa/dr < 0, and (Pipa/dr"^ < 0, and therefore ipa < 0, 
there can be no other zero for r < Vi. This proves that there is at most one zero for re (0, rj). 

But now we notice that near r = r,, equation (B9) has a singular solution (diverging as ln(ri — r)) and a unique 
regular solution of the form, 



V'(r) = ilj{r^) 



1 (r - n) + 0{r - nY 

niro-Ti) 



(BIO) 



Since Va is regular, it has the form (BIO), and this implies that close to r = the regular solution and its first 
derivative have opposite signs. This contradicts the result obtained under the assumption that there is a zero in 

< r < Ti. We conclude that (,/;„ docs not vanish for r € (0,ri). 

Similarly, we find that near r = ro, equation (B9) has a singular solution (diverging as ln(ro — r)) and a unique non 
vanishing regular solution. This implies that Ve{r) cannot vanish for r = Vg- 

Wc note in passing that for the extreme case Q = M, where = Tq, we have the exact (regular at the horizon) 

solutions, 

Mr) = ^(^r + 'M)ir-Mr^^^ ^^^^^ 

where C is a constant. 

Going back to (Bl), for any fixed £, given the solution (B6), a linearly independent solution is given by, 

tjjJr) = C-Veir) [ ^dy (B12) 

where C is a constant. It is easy to check that ^6 is regular and non vanishing in < r < r^, and, 

Mr) ~ Cir^ ; r^0+ 

Mr) ~ CiHn-r) ; r ^ (B13) 

where Ci and C2 are constants. We may obtain and expansion of this solution in powers of r using (B12), or directly 
from (Bl), 

4 , (16M-/3i)r5 5((£2^£-8)Q2+4M(12M-/3i))r6 
A{r) = r* + ^ „ h 



6Q2 42Q4 

{{ijt + 1){3()AI - rh) + - :m)M)Q'- + 3{)M-(32M - :iJi)')r' 

84Q6 



(B14) 



where we have set an arbitrary multiplicative constant so that the coefficient of is equal to one. 

For £ > 2 the integrals in (B12) cannot be computed directly, because that would require explicit expressions for 
the zeros of polynomials of degree larger that 4. We may, nevertheless, infer their general form as follows. We first 
notice that Ve{r) may be written in the form, 

Uk=ii-rk) 

where r j. arc the zeros of Ve{r), which, as indicated are simple. Therefore, since y^ — 2My + = {y — ri){y — rg), 
we should have, 

/ ^dy = Aln{ri-r)+Bln{ro-r)+C 

Jo {y^-2My + Q^){Ve{y))^ 

+ J2aMr-r,) + j:-^ (B16) 

fc=i fc=i ' 

where A, B, C, Uh and bk are constants that depend on M, Q, and rfe. The last term in (B16) may be written in the 
form, 

V _ Udr) .^.7. 
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where Ue{r) is a polynomial of order ^ — 1, or lower. Replacing in (B12), 

Mr) = ^Ve{r){Aln{n-r) + B\u{ro-r)+C) + ^Ue{r) 

. 1+2 

+-Ve{r)Y,a,\n{r-rk) (B18) 

^ k=l 

But we notice that '06 (r) is a solution of (Bl), which can be singular only at the regular singular points r = 0, r = ri, 
and r = ro, and that the zeros rk do not coincide with these points. Therefore, we must have Ofc = for all k, and the 
last term in (B18) vanishes identically. This result implies that, for any £, wc may construct algebraically the solution 
(B12) as follows. We first compute the coefficients of Pe{r) as indicated above, and then replace in (B18) leaving A, 
B, C and the coefficients of Ue arbitrary. Next we replace in (Bl) and impose the condition that '0 is a solution of 
that equation, and that tp ^ near r = 0. It can be checked that this procedure determines all the coefficients, up 
to an arbitrary multiplicative constant, a simple example being (83) for £ = 2. Since by construction these solutions 
satisfy the appropriate boundary condition at r = 0, the construction of the corresponding scalar zero modes, and the 
associated intertwining potential, is now a simple algebraic procedure. The resulting expressions are, unfortunately, 
very long and rather difficult to analyze in detail. In particular, wc have not been able to show explicitly that for 
< rc < Ti the scalar zero modes that satisfy the required boundary coiiditioii at r = are non vanishing everywhere 
in the interval < r < r^, as required for the regularity of the intertwining potential. Wc remark, nevertheless, that 
this appears to be the case in all the particular solutions analyzed numerically after assigning definite numerical values 
for the parameters, as in the example described in Section III - C. 
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